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, Abstract 

£^ . We calculate the (D(a cw ) electroweak radiative corrections to e + e~ — > tth at a electron-positron 

linear collider (LC) in the standard model. We analyze the dependence of the 0(a cw ) corrections 
on the Higgs boson mass rrih and colliding energy y/s, and find that the corrections significantly 
decrease or increase the Born cross section depending on the colliding energy. The numerical 
results show that the O(o; cw ) relative correction is strongly related to the Higgs boson mass when 
& \/s = 500 GeV, and for rrih = 150 GeV the relative correction ranges from —31.3% to 2.3% as the 

increment of the colliding energy from 500 GeV to 2 TeV. 
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I Introduction 

The main goal of the most experimental programs at present and future high-energy colliders is 
to search for Higgs boson, which is believed to be responsible for the breaking of the electroweak 
symmetry and the generation of masses for the fundamental particles in the standard model (SM) 
^ [2| . The fundamental particles acquire masses via the interactions with the ground state Higgs field 
in the SM Higgs mechanism. However, until now the Higgs boson hasn't yet been directly explored 
experimentally, except that LEP2 experiments provided a lower bound of 114.4 GeV for the SM Higgs 
boson mass at the 95% confidence level jH]. 

As we know that the Higgs search strategies depend largely on the suspected value of Higgs mass. 
Actually, it is most difficult to probe Higgs boson experimentally in the intermediate mass region 
(rrih ~ 100 — 200 GeV). In this mass region, the production mechanism with Higgs boson radiated 
from either a gauge boson or a fermion, is an important Higgs boson discovery channel. At e + e~ linear 
colliders and hadron colliders, the Higgs boson is searched via Bjorken process e + e~ (pp,pp) — > ffh, 
an intermediate Higgs boson is produced associated with a // pair. The coupling strength of the 
fermion-Higgs Yukawa coupling / — / — h is proportional to the fermion mass, i.e., 9ffh = mf/v, 
where v = (V2G F )-^ 2 ~ 246 GeV is the vacuum expectation value of the Higgs boson. Then 
we can see that the top quark Yukawa coupling g t f h is the largest one among all the fermion-Higgs 
couplings, e.g., g^ h ~ 0.5 to be compared for example with g?r ~ 4 x 10~ 4 . Therefore, the Higgs 
boson production via the process e + e~ — > tth is strongly enhanced by the top quark Yukawa coupling, 
and it can also be a basic mechanism for measuring the top quark Yukawa coupling. 
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Recently, a lot of efforts have been invested in improving the precise QCD theoretical corrections 
to the processes pp/pp — > tih + X |SJ |5J [Tj |S] • At a LC the cross section for e + e~ — > tih is small, 
about 1 fb for yfs = 500 GeV and mh = 100 GeV [HI 1101 ITT] . But it has a distinctive experimental 
signature and can potentially be used to measure the top quark Yukawa coupling in the intermediate 
Higgs mass region at a LC with very high luminosity. In Ref. |10j . S. Dawson and L. Reina found that 
the NLO QCD corrections increase the Born cross section by a factor of roughly 1.5 for e + e~ — > tih 
at yfs = 500 GeV and = 100 GeV. But at = 1 TeV, the corrections decrease the Born cross 
section and are relative small. These works show that the evaluation of radiative corrections is a 
crucial task for all accurate measurements of this process. 

In this paper we present the calculations of the full 0{a ew ) electroweak corrections to e + e~ — > tih 
in the SM. In section 2, we present our calculations of the full 0(a ew ) electroweak radiative corrections. 
The numerical results and discussions are presented in section 3. Finally, a short summary is given. 

II Calculations 

In this paper we use the t'Hooft-Feynman gauge. In the calculations of loop diagrams we adopt 
the definitions of one- loop integral functions of Ref. ^3]. The Feynman diagrams and the relevant 
amplitudes are created by FeynArts 3 |16j automatically, and the Feynman amplitudes are subsequently 
reduced by FORM ^7j. The numerical calculations of integral functions are implemented by using 
Fortran programs, in which the 5-point loop integrals are evaluated by using the approach presented 
in Ref. Q3]. 

The process e + e~ — > tih at the lowest level occurs through the Feynman diagrams of Fig.l. There 
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Figure 1 The tree- level Feynman diagrams for the e + e — > tth process. 

are two kinds of Feynman diagrams to the e + e~ — ► tth process at the tree-level. The first kind includes 
diagrams with Higgsstrahlungs from top or anti-top final states and the top Yukawa coupling is thus 
involved. The second kind consists of the diagrams with a Higgs boson radiated from a Z°-exchange 
s-channel or via Z° — G° — /i-interaction, and is independent of the top Yukawa coupling. We have 
checked that our numerical Born cross section is in good agreement with that in Ref . |10| 1121 11 4 j . 
The 0(a cw ) (one- loop level) virtual corrections to the process 

e+Oi) + e -(p2) -» t(h) + i(k 2 ) + h(k 3 ) (2.1) 

can be expressed as 

^virtual = CTo<5 vir tual = / d$ 3 Vtle (MoM^utd) . ( 2 ' 2 ) 

2 NV^ spin 

where cjq and A^o are the Born cross section and amplitude for e + e~ — > tth, respectively, d$3 is 
the three-body phase space element, the bar over summation recalls averaging over initial spins, and 
-^virtual ls the amplitude of one-loop Feynman diagrams and the corresponding counter-terms. Due 
to the fact that the electron-Higgs(Goldstone) Yukawa coupling is proportional to the electron mass, 
we do not consider the contributions of the one-loop diagrams which involve e — e — h(G°) vertex to 



the amplitude .M virtual- Therefore, the 0(a cw ) virtual corrections involve 975 loop diagrams which 
can be classified into self-energy (376), vertex (425), box (145) and pentagon (29) diagrams. As a 
representative selection, the pentagon diagrams are given in Fig. 2. 

The self-energy and vertex one-loop diagrams contain both ultraviolet (UV) and infrared (IR) 
divergences, while all the box and pentagon diagrams are ultraviolet finite and only contain IR diver- 
gences. To regularize the UV divergences in loop integrals, we adopt the dimensional regularization 
scheme |18| in which the dimensions of spinor and spacetime manifolds are extended to D = 4 — 2e. In 
this paper, we adopt the complete on-mass-shell (COMS) renormalization scheme |19j . in which the 
electric charge of electron e and the physical masses mw, m z, m-h, m t, m e are chosen to be the renor- 
malized parameters. The definitions and the explicit expressions of these renormalization constants 
can be found in Ref. ^H]- As we expect, the UV divergence contributed by the one-loop diagrams 
can be cancelled by that contributed by the counterterms exactly. Therefore, we get a UV finite cross 
section including 0(a ew ) virtual radiative corrections. 

The IR divergence in the process e + e~ — * tih is originated from virtual photonic corrections. It 
can be exactly cancelled by including the real photonic bremsstrahlung corrections to this process in 
the soft photon limit. The real photon emission process 

e + (pi) + e-(pa) -> t(h) + t(k 2 ) + h(k 3 ) + 7 (fc), (2.3) 

where the real photon radiates from the initial electron (positron) and the final top(anti-top) quark, 
can have either soft or collinear nature. The collinear singularity is regularized by keeping electron 
mass. In order to isolate the soft photon emission singularity in the real photon emission process, 
we use the general phase-space-slicing method [201- The bremsstrahlung phase space is divided into 

5 




6 



singular and non-singular regions, and the cross section of the real photon emission process (|2,3|) is 
decomposed into soft and hard terms 

Creal = °"soft + Chard = 0"o(<Wt + ^hard)- (2-4) 

For soft photons, ko < AE, we neglect the momenta of the radiated photons everywhere but in the 
singular propagators. By using the soft photon approximation, we find the contribution of the soft 
photon emission process is ^JJ l2*T] 

QePl _ QeP2 _ Qtfa Q t k 2 2 , 2 g v 

Pi ■ k P2 ■ k k\ ■ k k 2 • k ' 

in which AE is the energy cutoff of the soft photon and ko < AE <C \fs, Q e = 1 and Qt = 2/3 are the 
electric charges of the positron and top quark, ko = \J\k\ 2 + /i 2 is the photon energy. The integral 
over the soft photon phase space have been implemented, therefore, we obtain the analytical result of 
the soft corrections to e + e~ — ► tih which can be found in Refs. |19j and |21j . We checked numerically 
the cancellation of IR divergencies and verified the contribution of these soft photonic bremsstrahlung 
corrections leads to a IR finite cross section which is independent of the infinitesimal photon mass 
The hard photon emission cross section Uhard, with the radiated photon energy being larger than AE, 
is calculated by using Monte Carlo method. 

Finally the UV and IR finite total cross section including the full 0{a ew ) electroweak corrections 
reads 

a = a Q + Aa = <r + o- virtua i + a Ical = <t (1 + 8), (2.6) 
where 8 = ^virtual + <Wt + ^hard is the full electroweak relative correction of the order 0(et ew ). 



a ew f d 6 k 
dcr soft = -da — / -^r- 

^ J\k\<AE 2k 
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Ill Numerical results and discussions 

In our numerical calculations, we set a ew (0) _1 = 137.0359895, ray/ = 80.423 GeV, mz = 

91.188 GeV, m e = 0.511 MeV, = 105.7 MeV, m T = 1.777 GeV, m u = 66 MeV, m c = 1.35 GeV, 
m>t = 174.3 GeV, = 66 MeV, m s = 150 MeV, and mj = 4.3 GeV. The renormalization scale is 
taken to be Q = 2mt + rrih- Here we use the effective values of the light quark masses (m u and md) 
which can reproduce the hadron contribution to the shift in the fine structure constant a ew (mz) j-3j . 



Figure 3 The relative corrections of the 0(a ew ) order contributions to the e + e~ — ► tih cross section. 
The relative corrections of 5 v irtuai+soft and 5 hard, and their sum are shown as the functions of the 

soft photon energy cutoff AE/Ef,. 
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In Fig.3 we show the dependence of the 0(a ew ) relative correction to e + e~ — ► tih on the soft cutoff 
AE/Eb, assuming nih = 115 GeV and y/s = 500 GeV. As shown in this figure, both <5 v irtuai+soft(= 
Virtual + £ so ft) an< i <^hard depend on the soft cutoff AE/Eb, but the full 0(a ew ) electroweak relative 
correction 5 is cutoff independent within the range of statistical errors as expected. In the following 
calculations, the soft cutoff AE/Eb is fixed to be 10~ 3 . 



Figure 4 The Born and one-loop level cross sections for the process e + e — > tth as functions of the 

e + e~ colliding energy y/s. 



In Fig.4 we present the Born cross section ao and the one-loop level cross section a which include the 
full 0(a cw ) electroweak corrections, as functions of the e + e~ colliding energy y/s for the Higgs boson 



mass nth = 115 GeV, 150 GeV and 200 GeV, respectively. The colliding energy range (500 GeV ~ 
2 TeV) in the figure is accessible at future linear colliders, such as, TESLA [21] {^fs = 500 GeV), NLC 
M(y/s = 500 GeV), JLC (Vs = 500 GeV,) and CERN CLIC [ZZ] (1 TeV < y/s < 5 TeV). From 
this figure we can see that both curves of <ro and a for = 115 GeV reach maximal values at the 
position of y/s ~ 750 GeV, where the correction Act can reach 0.09 fb. For m;, = 150 GeV and 200 
GeV, both uq and a have their maximal values at y/s ~ 850 GeV and yfs ~ 950 GeV, respectively. 
The correction Act decreases as the increment of the Higgs boson mass m^. 



Figure 5 The dependence of the 0{a ew ) relative correction to e + e — ► tth on the e + e colliding 

energy y^i. 
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The dependence of the 0(a ew ) relative correction to e + e~ — > tih on the colliding energy y/s of a LC 
is displayed in Fig.5. For rrih = 115 GeV, the C(a ew ) corrections suppress the Born cross section in the 
energy region of y/s < 600 GeV, while enhance the Born cross section in the region of y/s > 600 GeV. 
The relative corrections can reach about —5.8%, —31.3% and —8.3% at y/s = 500, 500 and 600 GeV 
for rrih = 115, 150 and 200 GeV, respectively. For rrih = 150 GeV, the relative correction ranges from 
-31.3% to 2.3% as y/s running from 500 GeV to 2 TeV. The large correction for y/s = 500 GeV and 
rrih ~ 150 GeV shown in this figure comes from a threshold effect which diverges at the threshold. 



Figure 6 The Born and onedoop level cross sections for the process e + e — > tth as functions of the 

Higgs boson mass nih- 
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In Fig.6 we depict the Born cross section <ro and the one-loop level cross section a as functions of 
the Higgs boson mass rrih- As shown in this figure, Both the cross sections <7o and a decrease with 
the increment of the Higgs boson mass, and the cross sections <ro, <t and the one-loop level correction 
Act = a — ctq at yfs = 800 GeV are larger than those at y/s = 500 GeV and y/s = 2 TeV. 



Figure 7 The 0(a ew ) relative correction to the process e + e — > tth as a function of the Higgs boson 

mass m/i- 

In Fig.7 we plot the (D(a ew ) relative correction 5 to e + e~ — > tt/i as a function of m^. For y/s = 
500 GeV, the relative correction decreases from —5.8% at nih = 115 GeV to —31.3% at = 150 GeV. 
For y/s = 800 GeV and 2 TeV, the relative corrections decrease from 4.4% and 1.5% to —0.8% and 
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—4.7% as the increment of rrih from 115 GeV to 200 GeV, respectively. The 0(a ew ) electroweak relative 
corrections are not very sensitive to the Higgs boson mass in the range of 115 GeV < rrih < 200 GeV 
when y/s = 800 GeV and 2 TeV, but strongly depend on the Higgs boson mass when y/s = 500 GeV. 
In both Fig. 6 and Fig.7, we can see that each of the two curves of the total cross sections including 
the 0(a ew ) corrections at y/s = 800 GeV and 2 TeV, has two spikes at the vicinities of rrih = 2mw 
and rrih = 2m z, due to the threshold effects. 

After submitting this manuscript, we acknowledged another two papers appeared on this sub- 
ject [281 129j . The representative comparison with the calculation of Ref. 29^ is shown in Table 2 of 
29 (There we use the current mass values for m u and m^.). It shows that most of our numerical 
results of one-loop electroweak corrected cross sections agree with theirs within estimated error, but 
there are some discrepancies at the energy y/s = 2 TeV. 

IV Summary 

In this paper we calculate the full electroweak one-loop level radiative corrections to the process 
e + e~ — > tth at a electron-positron LC in the standard model. We analyze the dependence of the 
electroweak radiative corrections on the Higgs boson mass and colliding energy y/s, and find that 
the corrections increase or decrease the Born cross section in the Higgs boson mass range 115 GeV < 
rrih < 200 GeV, depending on the colliding energy. The numerical results show that the 0(a ew ) 
electroweak relative corrections can reach —31.3%, 4.4% and —4.7% at y/s = 500 GeV, 800 GeV and 
2 TeV, respectively. We also find that the full electroweak relative correction of the order 0{a ew ) is 
strongly related to the Higgs boson mass when y/s = 500 GeV, and for rrih — 150 GeV the relative 
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correction ranges from —31.3% to 2.3% as the colliding energy increasing from 500 GeV to 2 TeV. 
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Figure Captions 

Figure 1 The tree-level Feynman diagrams for the e + e~ — > tih process. 
Figure 2 The pentagon diagrams for the e + e~ — > it/i process. 

Figure 3 The relative corrections of the 0(a ew ) order contributions to the e + e~ — > ti/i cross 
section. The relative corrections of 5 v i r t ua i+soft and (5/^^, and their sum are shown as a function of 
the soft photon energy cutoff AE/Ej,. 

Figure 4 The Born and one- loop level cross sections for the process e + e~ — > tt/i as functions of 
the e + e~ colliding energy y/s. 

Figure 5 The dependence of the 0(a cw ) relative correction to e + e~ — > tf/t on the e + e _ colliding 
energy y/s. 

Figure 6 The Born and one- loop level cross sections for the process e + e~ — > tih as functions of 
the Higgs boson mass m^. 

Figure 7 The 0(a ew ) relative correction to the process e + e~ — > tt/i as a function of the Higgs 
boson mass m^. 
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